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Proof of the Double Bubble Conjecture 

By Michael Hutchings, Frank Morgan, Manuel Ritore, and Antonio Ros* 



Abstract 

We prove that the standard double bubble provides the least-area way to 
enclose and separate two regions of prescribed volume in M 3 . 



\J [ 1. Introduction 

Archimedes and Zenodorus (see [K, p. 273]) claimed and Schwarz [S] 
proved that the round sphere is the least-perimeter way to enclose a given 
volume in M 3 . The Double Bubble Conjecture, long believed (see [P, pp. 300- 
301], [B, p. 120]) but only recently stated as a conjecture [Fl, §3], says that the 
familiar double soap bubble of Figure 1, consisting of two spherical caps sepa- 
rated by a spherical cap or a flat disc, meeting at 120 degree angles, provides 
f->» | the least-perimeter way to enclose and separate two given volumes. 

Theorem (see 7.1). In M 3 , the unique perimeter- minimizing double 

bubble enclosing and separating regions R\ and R2 of prescribed volumes v\ 

and V2 is a standard double bubble as in Figure 1, consisting of three spherical 

caps meeting along a common circle at 120- degree angles. (For equal volumes, 

T ^ • the middle cap is a flat disc.) 

The analogous result in R was proved by the 1990 Williams College 
"SMALL" undergraduate research Geometry Group [F2]. The case of equal 
volumes in M 3 was proved with the help of a computer in 1995 by Hass, Hutch- 
ings, and Schlafiy [HHS], [Hu], [HS2] (see [Ml], [HS1], [M2, Chapt. 13]). In 
this paper we give a complete, computer-free proof of the Double Bubble Con- 
jecture for arbitrary volumes in R 3 , using stability arguments, as announced 
in [HMRR]. 

Rcichardt, Heilmann, Lai and Spielman [RHLS] have generalized our re- 
sults to M 4 and certain higher dimensional cases (when at least one region 
is known to be connected). The 2000 edition of [M2] treats bubble clusters 
through these current results. 
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Figure 1. The standard double bubble provides the least-perimeter 
way to enclose and separate two prescribed volumes. Computer 
graphics copyright John M. Sullivan, University of Illinois (http:// 
www.math.uiuc.edu/~jms/Images/) 

Previous results (see [M2, Chapts. 13 and 14]). Our strategy for proving 
Theorem 7.1 is to assume that a given double bubble minimizes perimeter 
and to use this assumption to deduce that the double bubble is standard. This 
strategy is valid only if we know that a perimeter-minimizing double bubble ex- 
ists. F. Almgren [A, Thm. VI. 2] (see [M2, Chapt. 13]) proved the existence and 
almost-everywhere regularity of perimeter-minimizing bubble clusters enclos- 
ing k prescribed volumes in M n+1 , using geometric measure theory. J. Taylor 
[T] proved that minimizers in M 3 consist of smooth constant-mean-curvature 
surfaces meeting in threes at 120-degree angles along curves, which in turn 
could meet only in fours at isolated points. An argument suggested by White, 
which was written up by Foisy [Fl, Thm. 3.4] and Hutchings [Hu, Thm. 2.6], 
shows further that any perimeter-minimizing double bubble in IR n+1 (for n ^ 2) 
has rotational symmetry about some line. 

Unfortunately, the existence proofs depend on allowing the enclosed re- 
gions Ri and R2 to be disconnected. The complementary "exterior" region 
could also a priori be disconnected. (If one tries to require the regions to be 
connected, they might in principle disconnect in the minimizing limit, as thin 
connecting tubes shrink away.) Hutchings [Hu] partially dealt with this compli- 
cation, using concavity and decomposition arguments to show for a perimeter- 
minimizing double bubble that both regions have positive pressure (see 4.1) 
and hence that the exterior is connected. Moreover there is a Basic Estimate 
(see §6) which puts upper bounds on the numbers of components of R\ and 
i?2j depending on the dimension n and the volumes v\, V2- 
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For equal volumes in R 3 , the Basic Estimate implies that both enclosed 
regions are connected. It can then be shown that a nonstandard perimeter- 
minimizing double bubble would have to consist of two spherical caps with a 
toroidal band between them (Fig. 8). Any such bubble can be described by two 
parameters, and Hass and Schlafly [HS2] used a rigorous computer search of 
the parameter space to rule out all such possibilities in the equal volume case, 
thus proving the Double Bubble Conjecture for equal volumes in M 3 . Earlier 
computer experiments of Hutchings and Sullivan had suggested that in fact 
no such nonstandard double bubbles were stable, and we confirm that in this 
paper, without using a computer. 

Our proof. In the present paper we consider arbitrary volumes vi, Vi 
in R 3 . We give a short proof using the Hutchings Basic Estimate that the 
larger region is connected (Proposition 6.2), and we use a stability argument 
(Proposition 6.5) to show that the smaller region has at most two components, 
as in Figure 2. (That the smaller region has at most two components can also 
be deduced from the Hutchings Basic Estimate using careful computation; see 
[HLRS, Prop. 4.6], [M2, 14.11-14.13].) 




Figure 2. A nonstandard double bubble. One region has two com- 
ponents (a central bubble and a thin toroidal bubble); the second 
region is another toroidal bubble in between. Computer graphics 
copyright John M. Sullivan, University of Illinois (http://www.math. 
uiuc.edu/~jms/Images/) 

To prove that an area-minimizing double bubble E is standard, consider 
rotations about an axis orthogonal to the axis of symmetry. At certain places 
on E, the rotation vector field may be tangent to E; i.e., the corresponding 
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normal variation vector field doiiS may vanish. The axis can be chosen so 
that these places separate E into (at least) four pieces (Proposition 5.8). Some 
nontrivial combinations w of the restrictions of v to the four pieces vanish on 
one piece and respect the two volume constraints. By stability, w satisfies a 
nice differential equation, and hence vanishes on more parts of £, which must 
therefore be pieces of spheres (Proposition 5.2). It follows that S must be the 
standard double bubble. 

The foregoing argument in the proof of Proposition 5.2 was inspired by 
Courant's Nodal Domain Theorem [CH, p. 452], which says for example that 
the first eigenfunction is nonvanishing. Other applications of this principle to 
isoperimetric problems and to the study of volume-preserving stability have 
been given by Ritore and Ros [RR], by Ros and Vergasta [RV], by Ros and 
Souam [RS] and by Pedrosa and Ritore [PR]. 

Open questions. We conjecture that the standard double bubble in M n+1 
is the unique stable double bubble. Sullivan [SM, Prob. 2] has conjectured that 
the standard /c-bubble in M n+1 (A; < n + 2) is the unique minimizer enclosing k 
regions of prescribed volume. This remains open even for the triple bubble 
in M 2 , although Cox, Harrison, Hutchings, Kim, Light, Mauer and Tilton 
[CHK] have proved it minimizing in a category of bubbles with connected 
regions (which a priori in principle might bump up against each other). 

One can consider the Double Bubble Conjecture in hyperbolic space M n+1 
or in the round sphere § n+1 . The symmetry and concavity results still hold 
[Hu, 3.8-3.10]. The case of § 2 was proved by Masters [Ma]. The cases of H 2 
and equal volumes in H 3 and in §> 3 when the exterior is at least ten percent of 
§ 3 were proved by Cotton and Freeman [CF]. 

There is also the very physical question in R 3 of whether the standard 
double bubble is the unique stable double bubble with connected regions. By 
our Corollary 5.3, it would suffice to prove rotational symmetry. In M 2 , Morgan 
and Wichiramala [MW] have proved that the standard double bubble is the 
unique stable double bubble, except of course for two single bubbles. 

Contents. Section 2 gives the precise definition of double bubble and a 
proof that there is a unique standard double bubble enclosing two given vol- 
umes. Section 3 provides variational formulas for our stability arguments. 
Section 4 gives some preliminary results on the geometry of hypersurfaces of 
revolution with constant mean curvature ( "surfaces of Delaunay" ) . Section 5 
uses stability arguments to show that a perimeter-minimizing double bubble 
in M ra+1 must be standard if one enclosed region is connected and the other 
region has at most two components. Section 6 proves the requisite component 
bounds for perimeter-minimizing double bubbles in R 3 . This completes the 
proof of The Double Bubble Conjecture, as summarized in Section 7. 
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2. Double bubbles 

A double bubble in W l+1 is the union of the topological boundaries of 
two disjoint regions of prescribed volumes. A smooth double bubble E C 
M n+1 is a piecewise smooth oriented hypersurface consisting of three compact 
pieces Si,S2 and So (smooth up to the boundary), with a common (n — 1)- 
dimensional smooth boundary C such that Si + So (resp. S2 — So) encloses a 
region R± (resp. R2) of prescribed volume v\ (resp. V2). None of these objects 
is assumed to be connected. The unit normal vector field N along S will be 
always chosen according to the following criterion: N points into R\ along dR\ 
and points into R2 along S2. We denote by a and H the second fundamental 
form and the mean curvature of S. Note that these objects are not univalued 
along the singular set C but they depend on the sheet Sj we use to compute 
them. We will also use the notation iVj, U{ and Hi to indicate the restriction 
of N, a and H to S;, i = 0, 1, 2. 

Since by Theorem 4.1 perimeter-minimizing double bubbles are smooth 
double bubbles (geometric measure theory automatically ignores negligible hair 
and dirt), throughout the rest of this paper by "double bubble" we will mean 
"smooth double bubble." 

A standard double bubble in R n+1 consists of two exterior spherical pieces 
and a separating surface (which is either spherical or planar) meeting in an 
equiangular way along a given (n — l)-dimensional sphere C. 

Proposition 2.1. There is a unique standard double bubble (up to 
rigid motions) for given volumes in M n+1 . The mean curvatures satisfy Hq = 
H\ — H.2- 

Proof. Consider a unit sphere through the origin and a congruent or 
smaller sphere intersecting it at the origin (and elsewhere) at 120 degrees as in 
Figure 3. There is a unique completion to a standard double bubble. Varying 
the size of the smaller sphere yields all volume ratios precisely once. Scaling 
yields all pairs of volumes precisely once. 

The condition on the curvatures follows by plane geometry for R 2 and 
hence for R n+1 (see [M2, Prop. 14.1]). □ 

Remark 2.2. Montesinos [Mon] (see [SM, Prob. 2]) has proved that there 
is a unique standard /c-bubble in M. n+1 for k < n + 2. 
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Figure 3. Varying the size of the smaller cap yields standard dou- 
ble bubbles of all volume ratios. Then scaling yields all pairs of 
volumes. 

3. Variation formulae 

In this section we will consider one-parameter variations { i Pt\\t\<e '■ ^ "" > 
W l+1 of a double bubble E C lR n+1 which are univalued along the singular 
set C and when restricted to each one of the pieces Ej are smooth (up to the 
boundary). Denote by X = dipt/dt the associated infinitesimal vector field at 
t = 0. Taking into account our choice of normal vectors to E it is a standard 
fact that the derivative of the volume of the regions R\ and R 2 are given by 

(3.1) -/" (JT.iVi)- [ (X,N ), and - /" (X,N 2 )+f (X,N ), 

Jj^i •-'So -'S2 •-'So 

respectively. On the other hand the first derivative of area of the bubble is 
given by 



=0,1,2 ^ 



where div^. is the divergence in Ej of a vector field in M n+1 . If {e.,} is an 
orthonormal basis of TEj and X is a vector field in M n+1 then div^.V = 
J2j \D e X, ej), where D is the Levi-Civita connection in M n+1 . As dw^ i X = 
div£jX T — nHi (X, iVj), where X T is the projection of X to TSj, the Diver- 
gence Theorem then implies the following well-known result. 

Lemma 3.1 (First variation of area for double bubbles). Let X C W l+1 
be a double bubble consisting of smooth hypersurfaces So, Si, S2, meeting 
smoothly along an (n— 1)- dimensional submanifold C. Then the first derivative 
of the area along a deformation (^t(S) at t = 0, where (ft is a variation with 
associated vector field X, is given by 

(3.2) - £ I nHi (X, Ni) - f (X, is + v x + v 2 ) , 

i=0,l,2 Ei C 

where iVj are the normal vectors to the smooth parts Ej of S and Vi are the 
inner conormals to C inside Ej. 
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Suppose that a double bubble £ is stationary for any variation preserving 
the volume of the regions R\ and R 2 - By Lemma 3.1 this is equivalent to 

(i) the mean curvatures Hi are constant, with —Hi + H2 + Hq = 0, and 

(ii) vq + vi + v 2 = along C. 

The mean curvature Hi (resp. H2) is called the pressure of the region Ri 
(resp. -R2). From (i) above, we get that if Hq > 0, then Ri has larger pressure 
than i?2- 

The functions Ui = (X, iVj) are the normal components of the variational 
field X. If the variation preserves volumes, from (3.1) they satisfy 



(3.3) / ui + / u = 0, / ui - I u = 0, 

and, since at the points of the singular set we have — Ni + N2 + Nq = 0, we 
get that 

(3.4) —ui + u 2 + u = along C. 

Now we follow the arguments in [BCE, Lemma 2.2] to show that any 
volume preserving infinitesimal variation is integrable. 

Lemma 3.2. Let £ C R n+1 6e a stationary double bubble. Given smooth 
functions m : Sj — ► R stic/i £/ia£ (3.3) and (3.4) are satisfied, there is a variation 
{(ft} of £ which leaves constant the volume of the regions enclosed by <£>*(£) 
and suc/t i/iai i/ie normal components of the associated infinitesimal vector field 
X are the functions Ui, i = 0, 1, 2. 

Proof. The boundary condition (3.4) allows us to construct a smooth vec- 
tor field Z on C such that (Z, iVj) = Uj, which can be extended smoothly along 
each Sj so that (Z,Ni) = m. Let {ipt} be a one-parameter variation of £ 
associated to Z (we can take Vt = V' + tZ). We choose nonnegative smooth 
functions /j : Sj — ► R, /« / 0, with compact support inside intSj, extended 
by to S. For t, si, «2 £ R close to 0, we consider the three-parameter 
deformation 

4>t + sifiN t 1 + S2 f 2 N t 2 , 

where Nj is the normal vector to ^(Sj). Let «j(t, si,S2), z = 1 2, be the 
volume of the deformed region Ri . Then 

£^(0,0,0)= / fj, 

which equals if i ^ j and is positive if i = j. From conditions (3.3) 

£(0,0,0)=* 
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Applying the Implicit Function Theorem we find smooth functions s\(t), s 2 (t) 
with Sj(0) = such that the volume of the regions Ri is preserved along the 
deformation. Let X be the vector field associated to this deformation. Using 
that the variation is volume preserving we get that s^(0) = 0. Hence the normal 
components of X are the functions U{. □ 

Now we wish to compute the second derivative of area for a variation of 
a double bubble keeping constant the volume of the two enclosed regions. 

Proposition 3.3 (Second variation of area for stationary double bubbles). 
Let £ C M n+1 be a stationary double bubble, and let (ft be a one- parameter 
variation with associated vector field X which preserves the volumes of R\ 
and R 2 . Then the second derivative of the area of <^t(S) at t = is given by 

(3.5) - J^u(Au + \a\ 2 u) - ^ L, Ui }7h7: + qiUi < > 

i=0,l,2' "' 



S ; _m,JC idVi 



where u = (X,N), U{ = (X,Ni), A is the Laplacian o/E, |<r| 2 is the squared 
norm of the second fundamental form, v^ is the unit inner normal to C inside 
Ej, and the functions qi are given by q± = (kq — n 2 )/V^, qi = ( — K i — K o)/v3 
and qo = («i + k 2 )/v3 with Ki = o~i{vi, v{), i = 0, 1, 2. 

Proof. First we recall that the derivative of the mean curvature H is given 
by 

(3.6) n^-(0) = Au+\a\ 2 u. 

dt 

To obtain the second derivative of area we differentiate (3.2) with respect 
to t. The derivative of the integrals over Ej in (3.2) equals 



(3.7) 



■f u{Au+\a\ 2 u)- V nHi — (f 
js i ^ 1 2 at t=0 \Jz t 



(X,N^ 



J* E . {X,Ni)Y Since the 



Let us see that the last sum vanishes. Let a,- = A 

dt t=o 

variation </?t(£) preserves volume, we obtain from (3.1) that a 1 + ao = and 
c-2 — «o = 0- As —Hi + H 2 + Hq = we conclude 

fliai + i? 2 a2 + #o«o = ao (-#i + fl" 2 + #o) = 0, 

which shows that the latter sum in (3.7) vanishes as we claimed. 

It remains to treat the boundary term in (3.2). Since uq + v\ + v 2 = 
on C, differentiating with respect to t we have 



d_ 

It 



t=o 



J (X, UQ + UX + v 2 )\ = J (X, D x (v + vi + v 2 ] 
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Equation (3.5) is then obtained from Lemma 3.6 below. To compute Dxv% the 
vector Vi has been extended as v\ along the integral curves of X, so that v\ is 
the unit inner conormal to y>t{C) in <^(£j). □ 

Remark 3.4. For bubbles in M 3 , the second variation formula and proof 
admit isolated singularities, such as tetrahedral soap film singularities. For 
bubbles in IR n+1 , C need only be piecewise smooth, including pieces meeting 
along an (n — 2)-dimensional submanifold. In addition, the second variation is 
insensitive to any sets of 7i n ~ 2 measure (see [MR, Lemmas 3.1, 3.3]). 

In a smooth Riemannian ambient manifold M ra+1 , the second variation 
has an additional term 

Ric(N,N)u 2 



involving the Ricci curvature in the normal direction N (see [BP, §7]). 

Remark 3.5. By approximation, the second variation formula (3.5) holds 
in a distributional sense (see (3.12)) for X piecewise C 1 or in H 1 . 

Lemma 3.6. Under the hypotheses of Proposition 3.3 we have 
(3.8) D x {v + v 1 + u 2 ) = E {f^J + »«i} N i = ~ {^r(°W» + ^(W 



where #2 = ^-{y\i v 2) ■, #0 = ^-iyiz v o) are the angles determined by the sheets of 
¥>t(E) along the singular set. 

Proof. Let Y = X c be the orthogonal projection of X to the tangent 
bundle TC. For each i we have X = Y + {X, i/j) i/* + (X, NJ Ni. As Y is 
tangent to C we have {Y,Dx(J2i u i)} = 0- We also have (ui,Dx^i) = for 
each i since |z/j| = 1. Moreover 

(D x Ui,Ni) = - (vi,D x Ni) = ai{vi,X % ) + — ^, 

where X 1 is the orthogonal projection of X to TEj. Hence 

f du- 1 

D x {v + vi + v 2 ) = E 1 IT. + CTi ^' ^) ^ X ' "*) + ai ( Y > Vi ^ f ^ 



i 



<9f,; 



Observe that [DyVi^vi) = 0, and that X)j.Dyi/j = ^VdZi^j) vanishes 
since Y is tangent to the singular set, where J2% u i 1S identically zero. Of 
course this implies Y^ii^Wi) = (Dy{J2i v i)) = 0- Hence we see that 

= Y^DyVi = E {( D YVif + (D Y Vi,Vi) Vi + (D Y v l ,N i ) TV*} 

i i 

= J2 (°y^ N i) N i = E < Y > v i) N i- 
i i 

Therefore we get that D x (v + v x + v 2 ) = J2i {§^ + «• (X, v t )} N,,. 
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Taking into account that 

(3.9) Ul = -L(N -N 2 ), U2 = -L(- Nl -N ), v = -L( Nl +N 2 ), 

we obtain that 

(X,ui) = — (u -u 2 ), (X,u 2 ) = —={-ui -uq), (X,u ) = —(u 1 +u 2 ). 

As -Ni + N 2 + iV = we have 

V3 ^2 Ki (X, Vi) Ni = V3^2qiUiNi + {(k u 2 - k 2 u ) Ni + (kwo - «o«i) N 2 

i i 

+(k 2 ui - KiU 2 )N } . 

The summand between brackets is a vector whose coordinates coincide, up to 
sign, with the determinant of the matrix 

/ N{ Kl m 

N% K 2 u 2 

where Nj = (Nj, . . . , Nf). But this determinant vanishes since — N\ + N\ + 
iVg = 0, —u\ + u 2 + uq = 0, and —k\ + n 2 + kq = 0. This last equality holds 
because —Hi + H 2 + i^o = 0, and —<j\(Z, T) + <72(Z, T) + <To(Z, T) = for any 
vector Z and T in TC. Hence the first part of (3.8) follows. 

To prove the remaining part we write v 2 = R(9 2 )vi, v§ = R(Qq)v\, where 
R(0) is the rotation in the plane spanned by vi and Ni given, in this basis, by 

R(9) = 

We have 

D x (vo + vi + vi) = ( Id + R(9 2 ) + R(0o))^r + ^(-^2) + ^HV>, 

and (3.8) follows since the first summand vanishes at t = 0. □ 

Remark 3.7. For a variation such that the angles of the sheets are pre- 
served, we have Dx(vo + vi + v-l) = (since vq + vi + v 2 = for all £), so by 
(3.8), the boundary term in the second variation formula (3.5) vanishes. 

Consider a stationary bubble E. We say that a function u : \J Ej — ► K 
defined on the disjoint union of the Ej is admissible if the restrictions Ui to the 
smooth pieces Ej of E, lie in the Sobolev space i? 1 (of functions in 1? whose 
gradient is squared integrable) and satisfy the boundary condition 

(3.10) ii i = u 2 + uq along the singular set C. 
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The space of admissible functions satisfying the zero mean value conditions 



'iUiVi 
1=0,1,2' " 

dui 

+ QiUi I Vi } , 



(3.11) / m+ u = u 2 - n = 

JSi jSo ^S2 J ^o 

will be denoted by .F(£). From the results at the beginning of this section, 
we see that admissible functions correspond to deformations of £ and that 
J~(E) are the infinitesimal variations of those deformations which preserve the 
volume of the regions R± and R 2 ■ The bilinear form on the space of admissible 
functions for the second variation of the area (3.5) will be denoted by Q, and 
it is given by 

(3.12) Q(u,v) = f {(Vu,Vv) -\a\ 2 uv} - ^ f q. 

"^ S i=0,l,2" C 

where i/j is the inner normal to C inside Sj and a, are the functions defined in 
the statement of Proposition 3.3. We will say that a (smooth) double bubble 
S is stable if it is stationary and Q(u, u) ^ for any u € ^"(S). We shall say 
that it is unstable if it is not stable. By Lemma 3.2 a perimeter-minimizing 
double bubble is stable. 

Lemma 3.8. Let £ be a stable double bubble and u £ J~(E) suc/i that 
Q(u,u) = 0. XTien u is smooth on the interior o/£j, i = 0, 1, 2, and t/iere 
exist real numbers Xq, Ai, A2, u>i£/i Ai = Ao + A2, suc/i i/iat 

Atij + |<r| tij = Aj, on Sj. 

Proof. The stability of £ implies that Q(u + to, u + to) ^ for any v £ J 7 
and i 6 I, Therefore Q(u, v) = and so, taking arbitrary functions with 
mean zero and support inside the interior of £j we conclude that the displayed 
equation holds in a distributional sense. From elliptic regularity, u is smooth 
on the interior of £j. □ 

A smooth admissible function u is said to be a Jacobi function if it corre- 
sponds to an infinitesimal deformation of £ which preserves the mean curvature 
of the pieces £&, and the fact that these pieces meet in an equiangular way 
along its singular set. By formulae (3.6) and (3.8), we have that u is a Jacobi 
function if and only if 

( Au + \a\ 2 u = 0, on £, 

(3.13) { (d Ul \ du 2 du . , n 

- q h fl«i = t\ \- Q2U 2 = q h go^o, along 6. 

Any Killing vector field Y of M n+1 gives a Jacobi function on £, u = (Y, iV). 
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Lemma 3.9. Let S C T, be a subdomain with piecewise- smooth boundary 
and u a Jacobi function on E which vanishes on dS (in particular we assume 
that all the U{ vanish at dS DC). If w is defined by 

{u, on S, 
0, on £ — S, 

then w is an admissible function and Q(w,w) = 0. 

Let S' C S be a second subdomain, with the same properties of S, and 
w' its associated admissible function. If the interiors of S and S' are disjoint, 
then Q(w,w') = 0. 

Proof. A Jacobi function u satisfies Q(u, v) = for any admissible function 
v by (3.12) and (3.13). Therefore the equalities Q(w,w) = Q(u,w) = prove 
the first assertion. The second one is trivial. □ 

4. Area-minimizing double bubbles and Delaunay hypersurfaces 

As described in the Previous Results section of the introduction, F. Alm- 
gren [A, Thm. VI. 2] (see [M2, Chapt. 13]) proved the existence and almost- 
everywhere regularity of perimeter-minimizing bubble clusters enclosing k pre- 
scribed volumes in M™ +1 , using geometric measure theory. Using symmetry, 
concavity, and decomposition arguments, Hutchings analyzed the structure of 
minimizing double bubbles. 

Theorem 4.1. 

(a) (after White, [Fl, Thm. 3.4], [Hu, Thm. 2.6]). An area- minimizing double 
bubble in R™ +1 (for n ^ 2) is a hypersurface of revolution about some line 
L. 

(b) ([Hu, Cor. 3.3]). In an area-minimizing double bubble, both enclosed 
regions have positive pressure. 

(c) ([Hu, Thm. 5.1]). An area-minimizing double bubble is either the stan- 
dard double bubble or consists of a topological sphere with a finite tree 
of annular bands attached as in Figure 4. The two caps are pieces of 
spheres, and the root of the tree has just one branch. All pieces are smooth 
(Delaunay) hypersurfaces meeting in threes at 120-degree angles along 
round (n — l)-spheres. 

Let O a connected component of the regions i?i or R2 in a nonstandard 
minimizing double bubble. Then either the smooth pieces in the boundary of 
£1 are all annuli or d£l is the union of two spherical caps D\ and D2 and one 
annulus Mq. In the first case we shall refer to fi as a torus component, and in 
the latter one as the spherical component. 
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Figure 4. Generating curve of nonstandard area-minimizing double 
bubble in R n+1 . 

Now we review some facts about hyper surf aces of revolution with constant 
mean curvature in JR™ +1 , known as Delaunay hypersurfaces (see [D], [E], [Ke] 
on M 3 and [Hs] on IR n+1 ). Let £ C M. n+1 be a hypersurface and assume that 
£ is invariant under the action of the group 0(n) of isometries of M n+1 fixing 
the xi-axis. The hypersurface £ is generated by a curve T contained in the 
xiX2-plane. The coordinates x\, X2, will be denoted by x, y, respectively. We 
parametrize the curve T = (x, y) by arc-length s. If a is the angle between the 
tangent to T and the positive x-direction we shall always choose the normal 
vector field TV = (sin a, — cos a). Then we have 

Lemma 4.2. The generating curve T of an O(n)- invariant hypersurface 
£ C M ra+1 with mean curvature H with respect to the normal vector N = 
(sin a, — cos a) satisfies the following system of ordinary differential equations 

(4.1) 



/ 
X 


= 


cos a, 




1 

y 


= 


sin a, 




a' 


= 


—nH + (n - 


, cos a 



y 

Moreover, if H is constant then the above system has the first integral 
(4.2) y n ~ x cos a - Hy n = F. 

The constant F in (4.2) is called the force of the curve T. Existence of the 
first integral is standard in the Calculus of Variations (see [GH, §3.4] and the 
references therein). For constant mean curvature surfaces see [P, pp. 138-139], 
with earlier reference to Beer and [KKS, §3]. 

From Lemma 4.2 we can obtain the following known properties. 

Proposition 4.3. Any local solution of the system (4.1) is a part of 
a complete solution T, which generates a hypersurface £ with constant mean 
curvature of several possible types (see Figure 5). 
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(i) // FH > then T is a periodic graph over the x-axis. It generates a 
periodic embedded unduloid, or a cylinder. 

(ii) // FH < then T is a locally convex curve and S is a nodoid, which has 
self -intersections. 

(iii) If F = and H ^ then T, is a sphere. 

(iv) If H = and F ^ we obtain a catenary which generates an embedded 
catenoid £ with F > if the normal points down and F < if the normal 
points up. 



(v) 
(vi) IfT, touches the x-axis, then it must be a sphere or a hyperplane. 



If H = and F = then T is a straight line orthogonal to the x-axis 
which generates a hyperplane. 



vn 



The curve V is determined, up to translation along the x-axis, by the pair 
(H,F). 



The generating curves of nodoids and unduloids are called nodaries and 
undularies. Since we shall often identify the curves and the generated hyper- 
surfaces we shall refer to them as nodoids and unduloids. 



Figure 5. The generating curves for Delaunay hypersurfaces: un- 
duloid, cylinder, nodoid, sphere, catenoid, hyperplane. 

Remark 4.4. Henceforth we shall use the following properties of gener- 
ating curves of Delaunay hypersurfaces. 

(i) Unduloids and nodoids have positive mean curvature with respect to the 
normal which points downward at the maximum of the y-coordinate. 

(ii) The nodoid is convex in the sense that the normal vector rotates mono- 
tonically This follows from equations (4.1) and (4.2). 
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Lemma 4.5 (Force balancing [KKS, (3.9)]). Assume that three generating 
curves Ti, i = 0, 1, 2 of hyper surf aces with constant mean curvature Hi and 
forces Fi meet at some point. If —H\ + H 2 + Hq = and —N\ + N 2 + iVo = 
at this point, then 

(4.3) -Fi + F 2 + F = 0. 

The lemma follows directly from (4.2). 

Lemma 4.6. Let £ be a nonstandard minimizing double bubble in R n+1 , 
as in Figure 4, and let R\ be the region of larger or equal pressure. Assume 
that the spherical component £1 is contained in R±. Let T\ be the generating 
curve of Mq = So n dQ,. 

Then the force of T\ is positive and T\ is an unduloid or catenoid and in 
particular a graph. 

Proof. Let T , Tq be the left and right circles in dQ. Consider the em- 
bedded curve T determined by Tq, r and Ti. Let F 2 be the third generating 
curve meeting F n Ti and T3 the one meeting Tq PI Ti. 

If the force of T\ is negative then T± is a nodoid with positive mean 
curvature since Hq = Hi — H 2 ^ 0. The graph T is convex and, since T meets 
L orthogonally, its total curvature equals n. At each one of the vertices T^nTi, 
T^nTi the inner angle of F is exactly 7r/3. By force balancing 4.5, both T 2 and 
Fz have positive force and they are unduloids. Since R 2 has positive pressure 
both T 2 and T3 are inward graphs with respect to Ti (i.e., the exterior region 
lies above Ti and above T 2 ). Hence the two circular arcs T , Tq have angular 
measure larger than ir/3. So the total curvature of T is larger than 47r/3, which 
is a contradiction. 

If the force of Ti is then Ti is part of a circle or of a line orthogonal 
to the axis of revolution L. The former possibility is discarded by the same 
argument used for nodoids. The latter possibility is clearly not possible. 

Hence the force of Ti is positive and Ti is an unduloid or catenoid and 
graph. □ 

By a similar argument to the one used in Lemma 4.6 we obtain 

Lemma 4.7. Let T, be a double bubble of revolution such that both regions 
have positive pressure. Then it is not possible that E contains pieces of spheres 
Tq, Tq, Ti, T 2 , T3 as in Figure 6, with Ti C So- 

Lemma 4.8. Let £ be a nonstandard minimizing double bubble in M 3 . 
Let 9i be the subtending angle of the spherical caps Di as in Figure 7. 

(i) If 0\, 9 2 < 7r/6 then 6\ = 2 and Mq is symmetric with respect to a 
plane orthogonal to the axis of revolution. 

(ii) If6 1 < vr/6 < 6 2 < tt/3 then 6 2 >--6 1 . 

o 
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Figure 6. 




Figure 7. The spherical component £1 

Proof. Assume that the spherical component O is contained in R\. By 
scaling, normalize H\ = 1. 

Let T\ be the generating curve of Mq parametrized from the left to the 
right. As oii = 6i — tt/6 we get from (4.2) that the force Fq of T± is given by 
0(01 ) = 5(02), where 



9(0) 



1 \ f^x 

Hn ) sin 6 -\ cos 6 sin ( 

2 V 2 



Since H = 1 — H 2 < 1 by Proposition 4.1, we get g'(9) ^ -^ sin 2(9 + 
■^ cos 20 > 0, and so the function g is strictly increasing in [0, 7r/6]. Hence 6±, 
9 2 ^ 7r/6 implies 6\ = 9 2 . Moreover the endpoints of T\ have the same height 
Hi and the same angle aj. This proves (i) by the uniqueness for solutions to 
(4.1) with respect to the initial conditions. 

To see (ii) let 

h{6) = — sin 2 6> + — sine cos 9 = g{9) - (1 - H ) sin 2 9. 

It is easily proved that the function h is symmetric with respect to ir/6 
and increasing in [0, 7r/6]. Thus we have 



9 



TV 



1 -/To)si^(|-M+M|-0 2 



< (1 - H ) sin 2 9 2 + h{9 2 ) = g{9 2 ) = g(6 1 ) 
and, as g is increasing in [0, vr/6], we get (tt/3) — 9 2 < 0±, as we wished. □ 
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Remark 4.9. If £ were n-dimensional, then the force of T\ would be 
given by Fq = g{6) sin n ~ 2 6 and so Lemma 4.8 works in arbitrary dimension. 

5. Separation and instability 

Let £ C R n+1 a stationary double bubble of revolution whose axis L is 
the xi-axis with generating curve V C {(xi, x 2 )| x 2 ^ 0} consisting of circular 
arcs To meeting the axis and other arcs Tj meeting in threes, with interiors Tj 
(see Figure 9). The bubble £ is invariant under the action of the group 0(n) 
of orthogonal transformations in R n+1 which fix the xi-axis. We consider the 
map / : r — L — > L U {oo} which maps each p £ T — L to the point L(p) n L, 
where L(p) denotes the normal line to T at p. If L(p) does not meet L, we 
define the image of p as f(p) = oo. Note that / is multivalued at the endpoints 
of the arcs Tj, where three of them meet. We will use the notation iA and iB 
for the image under / of the endpoints of Tj. 

Remark 5.1. Using (4.1) and (4.2), we find that the derivative of / with 
respect to arc length is given by /' = n _f — 2—. In particular, / is increasing 
as we move to the right along an unduloid or the concave up portion of a 
nodoid, decreasing as we move to the right along the concave down portion of 
a nodoid, and of course constant on spheres and vertical hyperplanes. Hence 
/ is locally injective on any Delaunay curve with nonzero force. 

Proposition 5.2. Consider a stable double bubble of revolution E C 
R n+1 , n ^ 2, with axis L. Assume there is a finite number of points {p±, . . . ,p^} 
in [JiTi with x = f(pi) = • • • = f{pk) which separates T. Assume further that 
{pi, . . . ,pk} is a minimal set with this property. 

Then every connected component of So, Ei and £2, which contains one 
of the points pi is part of a sphere centered at x {if x G L) or a part of a 
hyperplane orthogonal to L (if x = 00). 

Proof. First suppose that x £ L and take, after translation, x = 0. The 
1-parameter group of rotations 

<Po(xi, . . . ,x n+ i) = (cos(9xi + sin#x 2 , -sin(9xi + cos<9x 2 ,x 3 , • • • ,x n +i), 
8 £M, determines a Jacobi function on the bubble, u : E — ► R, given by 

d 

where N(p) is the unit normal vector of E at p, {e±, . . . , e n+ i} is the standard 
orthonormal basis of M n+1 and det denotes the Euclidean volume element. We 
define here Mq = E n {x 2 = 0}. By the symmetry of E, if p € Mq, then the 
vector N(p) also lies in the hyperplane x 2 = and therefore u vanishes on Mq. 



u{p) ^de 



<Pe(p),N(p)\ = -det(p,N(p),e 3 , . . . ,e n+1 ), 

1 / 
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On the other hand, if we take p in / -1 {0}, then the vectors N(p) and p 
are collinear. Using again the invariance of S with respect to 0(n), we get 
that u vanishes on the orbit M{p) of p under the action of 0{n) (note that 
M(p) is a hypersurface of £). 

As the points pi, . . . ,pk separate the curve T, the set M(pi)U. . .UM^JU 
Mq is a hypersurface of the bubble contained in u _1 {0} which separates X in 
at least four connected components. In fact, as the set {pi, . . . ,pk} is minimal 
among the subsets of / _1 {0} satisfying the separation property, it follows that 
£ — [M(pi) U . . . , M(pk) U Mo] has exactly four components Ai, . . . , A4 and 
that each one of the sets M(p±), . . . ,M(pf c ),Mo meets the boundary of each 
one of these four components. 

Consider the functions yW, i = 1, . . . , 4, on £ given by 



v 0) 



u, on Aj, 

0, on £ — Aj. 



Then v^ % > are admissible and we can find scalars a±, 02, 03, not all equal to 
zero, such that v = J2i=i a i v verifies the mean value conditions (3.11), so 
that v G -^(X). By Lemma 3.9, 

Q(v,v) = Y,alQ(v {i) ,v {i) ) = 0. 

Since u is a Jacobi function, 

(5.1) Av + H 2 t> = 

on S \ [M(p\ ) U • • • U M{pk) U Mo] . By our stability hypothesis and Lemma 3.8, 
equation (5.1) holds on all of S. 

Fix i and let S be the connected component of a smooth piece of E which 
contains the point pi. As pi lies in the interior of S, the four domains Aj 
meet the interior of S. As v vanishes onSfl A4, from the unique continuation 
property, we conclude that v = on S. Hence u = on S n Aj, for any 
j £ {1, 2, 3} such that aj / 0. As such j exists we conclude that u = on S 
again from the unique continuation property. Thus the 1-parameter group of 
rotations ifg preserves S. Since S is rotationally symmetric around the Xi-axis, 
we conclude that this component is a part of a sphere centered at the origin. 

This finishes the proof of the proposition if a; is a finite point of the axis L. 

It remains to prove the result when x = 00. In order to prove it we re- 
peat the argument by considering, instead of the rotations </?#, the 1-parameter 
group of translations Tq(x\, . . . , x n+ i) = {x\, X2 + 0, . . . , x n+ \) and its associ- 
ated Jacobi function u(p) = (N(p),e2)- □ 

Corollary 5.3. There is no stable double bubble of revolution in R n+1 
in which the graph structure is the one in Figure 8. 
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Figure 8. There is no stable nonstandard double bubble with con- 
nected regions. 

Proof. Assume first that the line equidistant from the two vertices inter- 
sects the axis L in a point p. Then Ti and T2 each has an interior point farthest 
from or closest to p, so that p £ /(Ti) n /C^). By Proposition 5.2, Ti and T2 
are both spherical (centered on the axis), which is impossible. 

If the equidistant line is horizontal T±, T2 each has an interior point far- 
thest left or right, so that 00 £ /(Ii) n /O^). By Proposition 5.2, Ti and T2 
are both vertical, which is impossible. □ 

Remark 5.4. When n = 2 and the volumes are equal, Hutchings [Hu, 
Thm. 5.1, Cor. 4.4] showed, as described in our Section 6, that any nonstandard 
minimizing bubble satisfies the hypotheses of Corollary 5.3. Therefore in this 
case the minimizing solution is the standard bubble. This fact was first proved 
by computer analysis by Hass and Schlafly [HS2]. 

Corollary 5.5. Consider a stable double bubble of revolution in M n+1 
in which both regions have positive pressure. Assume that one of the regions 
i?2 is connected, that the other one R± has two components and that the graph 
structure is the one in Figure 9. 

Then there is no x G L such that f -1 (x) — To contains points in the 
interiors of distinct Tj which separate T. 




Figure 9. A candidate double bubble of three components. 

Proof. There must be points on T±, T2, or T3. By Proposition 5.2, one 
of them is spherical. By force balancing 4.5, all three are spherical, which is 
impossible by Lemma 4.7. □ 

We give the following version of Corollary 5.5, although we will not use it 
here. 

Corollary 5.6. For a stable double bubble of revolution, if f is not 
injective on the interior ofTi, then Ti is a circular arc or vertical line. 
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Proposition 5.7. Consider a minimizing nonstandard double bubble in 
R n+1 , which is necessarily rotationally symmetric around an axis L. 

Then there exists no x G L such that / -1 (x) — Tq contains points in the 
interiors of distinct Tj which separate T. 

Proof. Since every component borders the exterior, a separating set must 
cut the outer boundary of some component. By force balancing 4.5, every 
piece of the outer boundary of this component is spherical. If f~ 1 (x) cuts two 
pieces of the outer boundary, then these are pieces of spheres with the same 
center and the same mean curvature, and hence the same sphere. The portion 
of the bubble between these two pieces can then be rolled around the sphere, 
without changing perimeter or enclosed volume, until it touches some other 
part of the bubble, resulting in a bubble which is not regular, and hence not 
minimizing. So it cuts an inner boundary (part of So). By force balancing 4.5, 
each end of the inner boundary meets two other spheres, which contradicts 
Lemma 4.7. □ 

Proposition 5.8. There is no stable double bubble of revolution in W l+1 
in which both regions have positive pressure, the region of smaller or equal 
pressure Ri is connected, the other region R\ has two components, and the 
graph structure is the one in Figure 9. 

Proof. Suppose there were. To are spherical. T\ is an unduloid or catenoid 
and graph by Lemma 4.6. By force balancing 4.5, T2 and T3 are (convex) 
nodoids. Since the top, third component has larger pressure, T4 must be a 
(convex) nodoid, catenoid, or vertical line, unless it is upside down (which 
cannot occur in the principal cases of Figure 14). (Here by "convex" we just 
mean that the tangent vector rotates monotonically.) 

We focus on the third component and its two vertices ^245 and V345. For 
the simplest case when all the curves are graphs as in Figure 10A, then the 
images iA and iB under / of the left and right endpoints of Tj satisfy 
AA < 2B < 5A and 5B < 3A < AB. 

This remains true as a vertex rotates until one of the three tangent vectors goes 
vertical. (The borderline position with 5A = 00 may be considered an extreme 
position of either case; in the proof we consider it part of the second case 
eliminated.) Rotating ^245 counterclockwise one notch as in Figure 10B yields 
instead 5 A < A A < 2B. To avoid giving T4 or Tr, two vertical tangents contrary 
to Corollary 5.6, the two vertices must be rotated in the same direction, say 
counterclockwise. Suppose that ^245 is rotated mi notches counterclockwise 
and that W345 is rotated mi notches counterclockwise. Then rri\ < 2, or T2 
(where R2 is a convex region by positive pressure) could not meet the circle 
To at 120 degrees (see Figure 11). Also mi < 3, or T4 would go vertical twice 
(see Figure 12), contrary to Corollary 5.6. 
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A 



B 



Figure 10. Rotating vertex U245 one "notch" (counterclockwise) 
means turning one tangent (here T§) past vertical. (Rotating an- 
other notch would turn T2 past vertical.) 




Figure 11. If V245 turns three notches, T2 cannot be convex. 




Figure 12. If V345 turns four notches, T4 goes vertical twice. 

Next consider cases with 771,2 = 2 or 777,2 = 3, as in Figure 13. T3 is not a 
graph, /(Ts) = [00, 3^4) U (3B, 00], and by stability Proposition 5.2 gives that, 
for T3, 3 A < 3B. We then have 3B < IB, or else T3 would go vertical a second 
time near 3B, contradicting Corollary 5.6. Therefore IB is contained in /(T^), 
which contradicts Corollary 5.5 for T\ and T3. 




Figure 13. If 7772 = 2 or 7772 
/(r 3 ) = [oo J 3A)U(3B,oo]. 



3, then T3 is not a graph and 
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Now remain only the cases < m\ ^ 2, < m-i < 1 of Figure 14. (It is 
easy to see that T4 cannot be a vertical line.) We claim that 



(5.2) 



3A < /(r x ) 



This is easy if T3 is a graph, since consideration of i>i3 shows that 2>B < IB, 
and then f(T^) < f(Ti) by Corollary 5.5 for T\ and T3. V13 can rotate only 
clockwise one notch to keep the stem part of a circle and T± a graph. Now T3 
is not a graph and /(T^) includes [oo,3j4). By Corollary 5.5 for Ti and T3, 
3A < f(Ti) as claimed. 




m 2 



m 1 



Figure 14. The six principal cases to be eliminated 

For the cases (0,0), (0,1), a similar argument shows that f(Ti) < 2B. 
Consideration of the vertices v 245 and V345 leads to the conclusion 



(5.3) 



5B <3A< /(fi) <2B < 5A 



Since the net angle #5 through which T$ turns satisfies #5 < 180 degrees, 
obviously 4A < 3A. For the case (0,1), where /(r4) contains (4.A, 00], this 
puts 3^4 in f(T4) n fi^s), a contradiction of Corollary 5.5 for T3, T4, T^. For 
the case (0, 0), consideration of ^345 shows that 3^4 < 4B and leads to the same 
contradiction. 
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Next we consider the cases (1,0), (2, 0). Since 5B < 3A, 3A is contained in 
/(r5). Since 3A < 4B, by Corollary 5.5 for l - ^, I^, T$, we must have 3^4 < 4A 
In particular, 6§ > 180 degrees. Of course by Corollary 5.6 for T§, 5A < 5B. 
Moreover T§ leaves ^345 above the horizontal. Now Corollary 5.10 implies that 
3 A > 4A, a contradiction. 

Similarly for the final cases (1, 1) and (2, 1), 3A is contained in f(T$). Since 
/(r4) includes (4A, 00], by Corollary 5.5 for 1^3, I^, T§, we must have 3^4 < 4A, 
an immediate contradiction in case (1, 1). In particular, #5 > 180 degrees, and 
5A ^ 5-B. If T§ leaves 1*345 at or above the horizontal, Corollary 5.10 yields 
the contradiction 3A > AA. If on the other hand T§ leaves ^345 below the 
horizontal, then the downward normal n to I^ at V345 is counterclockwise from 
the downward tangent to 1^ at W12 (and hence from every downward tangent 
to I"2) and hence counterclockwise from the downward normal to T\ at v±2- 
Since I\ is convex, n stays to the right of I^ and 1A < 3A, a contradiction of 
(5.2). □ 

Lemma 5.9. Given points A and B, consider two points D, E on the 
same side of AB subtending the same angle 9 as in Figure 15. Then ZCDE = 
ZABC. 



Proof. Since ZBCA = ZDCE and ACE ~ BCD, ABC ~ CDE. 



a 




Figure 15. ZCDE = ZABC 

Corollary 5.10. In cases (1,0), (2,0), and (2,1) of Figure 14, suppose 
that the net angle #5 through which T§ turns exceeds 180 degrees, that T^ leaves 
t>345 at or above the horizontal, and that 5 A < 5B. Then 4A < 3 A. 

Proof. Let = 6$ — 180 > 0. Apply Lemma 5.9 with A = V345, B = V245, 
AD _L T§ and AE _L T3 at W345, BD _L T§ and BE _L T4 at ^245; then 
ZADB = ZAEB = 9. Since T4 is strictly convex (it cannot be a vertical line 
because 5,4 s^ 5B), ZABC > 90. By Lemma 5.9, ZCDE = ZABC > 90. Since 
by hypothesis T$ leaves U345 at or above the horizontal, DE heads downward. 
(In these cases, Figure 15, in which AD is vertical, is rotated clockwise by an 
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amount less than 90 degrees, strictly less by Corollary 5.6 since #5 > 180.) 
Since 5A ^ 5B, D lies on or below the horizontal axis. Hence E lies below the 
horizontal axis and AA < 3 A. □ 



6. Estimates on the number of components 

In this section we prove that for a minimizing double bubble in R 3 , the 
larger region is connected, and the smaller region has at most two components. 

We begin by recalling a version of the Hutchings Basic Estimate. Let 
A(v) denote the volume of a sphere in M 3 enclosing volume v, and let A(v±, v 2 ) 
denote the area of the standard double bubble enclosing volumes v\ and v 2 in 
IR 3 . We then have: 

Proposition 6.1 ([Hu, Thm. 4.2]). In an area- minimizing double bubble 
enclosing volumes v\ and V2 in M 3 , suppose that R 2 contains a component of 
volume Xv 2 . Then 

2A(v!,v 2 ) ^ \- 1/3 A(v 2 ) + A(v 1 )+A(v 1 +v 2 ). 

This inequality places a lower bound on A in terms of v\ and v 2 . Hutchings 
[Hu, 4.4, 4.5] calculated that when v\ = v 2 , the lower bound is greater than 
1/2, so both regions are connected; and when v\ » v 2 , the lower bound 
approaches 2/5, so the smaller region has at most 2 components. More recently, 
Heilmann et al. [HLRS, Fig. 8] ([M2, Fig. 14.11.1]) used a computer to plot 
the lower bound on A as a function of the ratio v 2 /vi, and found that it 
is apparently increasing. This would imply that the larger region is always 
connected and the smaller region always has at most two components. However 
this observation has not been rigorously proved, because the function A(v 1, v 2 ) 
is difficult to work with. Thus we will use different methods to prove the above 
connectedness results. 

Proposition 6.2. In a minimizing double bubble in M 3 , the region with 
larger or equal volume is connected. 

Proof. By rescaling, we may assume that the two volumes are 1 — w and w. 
Hutchings [Hu, Thm. 3.5] showed that if w < 1/3, then the larger region is 
connected. (This is true in higher dimensions as well.) 

For w > 1/3, to prove that the region of volume 1 — w is connected, it is 
enough to show that the inequality of Proposition 6.1 fails for A = 1/2; i.e. 

2A(w, 1 - w) < 2 1/3 A(1 -w) + A(w) + ,4(1). 

We observe that by Lemma 3.1, dA(w, 1 — w)/dw > for w < 1/2, be- 
cause we can continuously deform one standard double bubble to another, 
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and the smaller region has larger pressure. Thus A(w, 1 — w) < ,4(1/2, 1/2). 
It is straightforward to compute that ,4(1/2,1/2) = 2" 4 / 3 3A(l), and A(v) = 
v 2 ' 3 A(l). Thus it will suffice to show that 

2-V3 3< 2V3( 1 _ w )2/3 +w 2/3 + L 

Since this holds at w = 0.10 and w = 0.63 and the right-hand side is concave, 
it holds for 0.10 ^ w < 0.63. In particular, it holds for 1/3 ^ w ^ 1/2, as 
desired. □ 

Remark 6.3. An alternative proof is given by Heilmann et al. [HLRS, 
Prop. 2.5]. Actually as in the proof above, any region with at least 37% of the 
total volume is connected. 

Lemma 6.4. In a minimizing double bubble in IR n+1 enclosing two un- 
equal volumes, the smaller region has larger pressure. 

Proof. Consider the function A{v , 1 — v) giving the least area enclosing 
and separating regions of volume v, 1 — v. By Hutchings [Hu, Thm. 3.2], 
A is strictly concave and of course symmetric about v = 1/2. Moving the 
separating surface (of mean curvature Hq) of an area-minimizing bubble we 
have dA/dv = uHq, and the left and right derivatives of A must satisfy 

A' R < nH s; A' L . 

Consequently Hq is positive for v < 1/2 and negative for v > 1/2. In other 
words, the smaller region has larger pressure. □ 

The following Proposition 6.5 shows that small region with three compo- 
nents is unstable, by using a volume-preserving linear combination of (non- 
volume-preserving) variations of each of the components. 

Proposition 6.5. In a minimizing double bubble in M 3 , the region with 
smaller or equal volume has at most two connected components. 

Proof. Assume that the volume of R± is less than or equal to the volume 
of i?2- By Proposition 4.1 we obtain Hi, Hi > 0. By Lemma 6.4 we get 
H ^0. 

Recall that Ki = o~i{i>i,Vj), and let Cj = o~i(T,T), where T is the unit 
tangent vector to the singular curve C. So 2Hi = m + Cj. 

We consider an admissible function u invariant by the one-parameter 
group of rotations of S. The functions Ui are locally constant over C. If 
we apply (3.12) to u, adding and subtracting CiU? in the boundary term, we 



484 M. HUTCHINGS, F. MORGAN, M. RITORE, AND A. ROS 

see that the second variation form satisfies 

Q(u,u) = W {|Vu,| 2 - \<Ti\ 2 « 2 } 
i ^ s * 

- A |^(F - H 2 ) u\ + (-ili - Ho) u\ + (ili + ff 2 ) ^o 

+ -/5 / ( C - c 2) Ml + (~Cl - C ) U 2 2 + (Ci + C 2 ) Uo- 

Taking the scalar product with L>tT in the formulae (3.9) we have 

4 = ~^( c o- c 2)> «p = 77f(~ Cl ~ c °)' K ° = 77|( c i+ c 2) 5 

where k* stands for the geodesic curvature of C inside E$ (with respect to the 
conormal z/j). So we have 

(6.1) Q(u,u) = E/ {|V^| 2 -h| 2 uf} 

- -^ ^(tf - #2) w? + (-#i - H ) u 2 2 + (Ifi + H 2 ) ul 

Consider a connected component f2 of R\. Let Mj = S^ndO, and let (7* = 
C n dfi. We want to find an admissible function u such that Q(u, u) < with 
support inside dQ. Then if R\ had three connected components, some nonzero 
linear combination would preserve the volumes and yield a contradiction. 

We define the function 

f 1, on M UMi, 
I 0, elsewhere in S. 



Then (6.1) gives 

2 

/ i^r - ■ 

i=0,l ' 



(6.2) Q(v,«) = -E / l^| 2 -^f/ (Ho + H^+f K + 4). 



Since |<7j| 2 = 4H?—2Ki, from (6.2) and Gauss-Bonnet J M . Kj = 2nx(Mi) 
JdM t K l we obtain 



(6.3) Q(v,v) 



J2 ^ X (Mi)-J M 4H?y^= JjHo + m)- Jj K g + K l). 



Assume first that f2 is a torus component, so that its boundary is a union 
of annuli. Adding (6.2) and (6.3) and taking into account that x(-^i) = 0) we 
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eliminate the geodesic curvature to obtain 

2Q(v,v) = -J2 I {N 2 + 4i^}-4^/ (flo + fli)<0, 

as desired. 

We now assume that O is the spherical component, so that M\ is the 
union of two spherical caps D±, D2 and an annulus Mq, as in Figure 7. As 
Mq is a graph by Lemma 4.6 we conclude < 8i ^ -^-, where 6i is the angle 
determined by D^. By scaling we may assume that the spherical caps have 
mean curvature H\ = 1. Using Gauss-Bonnet we get that 

(6.4) A(M 1 )= f K 1 = A7r- f k\. 

J Mi JC* 

Since vq = (—1/2) ^1 + (\/3/2)JVi, taking the scalar product with DtT we 
have 

(6.5) ,° = -^ + # 

From (6.3), (6.4) and (6.5) we obtain, taking into account that %(Mi) = 2 
and discarding the summands containing Hq, 

Q( ViV ) <en+ 'LA{M 1 )-^=L(C*). 

As A(Di) = 2ir (1 — cos#j) and L(dDi) = 2irsm9i, we have 

(6.6) Q(w,v)<27r{-4 + /i(0i) + &(0 2 )}, 

where /i(#) = | ( cos# -U sin# ), which is decreasing in [0, 2ir/3]. Thus if 6\ 

or #2 is greater than or equal to tt/2, we have 

Q(v,v) ^ 2ir{h(0) + h(n/2)} < 0. 

Assume now that both 8±, 62 < 7r/2, and consider the function 

cos# 

7-, m A, 

_ COS 0j 

"' ( 1, m M , 

0, elsewhere in S. 

As v and u; differ only on Mi we obtain from (3.12) 

(6.7) Q(w,w)= [ (\Vw\ 2 -2w 2 ) +2 / 1 + Q(w,u). 

■/Mi v y J M\ 

By direct computation we get 

Vw\ 2 - 2w 2 ) = -2tt Y ^|i, / l = 2vr V (1-cos^), 
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which combined with (6.6) and (6.7) yield 

(6.8) Q(w,w)^2ir{g(8 1 )+g(e 2 )}, 

where g is given by 

ra\ 3 a 7 ■ r sin ^ 
g{6) = - cos 9 — sin ( 



2 2\/3 cos 9 ' 

The function g is strictly decreasing in [0, 7r/2] since it is the sum of three 
decreasing functions. As g(0) = | > 0, g(ir/6) = 0, and g(ir/3) = —5/2 we 
conclude 

IT IT 

Q(w,w) < if either both 61,62 > — or some 9i ^ — . 

6 3 

We finally consider the remaining cases in which at least one of the angles 
6i is less than or equal to 7r/6 and both are less than tt/3. 

Case 1. 6»i < f < 6 2 < f . 

Observe that g is concave in [0, 7r/3] since 

7 / 1 „\ / sin 2 6* „ sin 3 ■ 



g"(9) = -- [cose -^= sm9) - 3 ^ + 2 — ^- < 0. 

2 V V3 J \ cosd cos 2 9) 

By Lemma 4.8 we know that ^ — 9\ < 9 2 - As g is decreasing and concave 

±-Q( w , w )^g(9 1 ) + g(9 2 )<g(9 l )+g(^-9 1 ^ <^g(^j =0. 

Case 2. 6>i, 2 < f . 

By Lemma 4.8 Mq is symmetric with respect to a plane orthogonal to the 
line of symmetry. So if k\ = o\{y\, v\) ^ we get from (3.12) 

Q(v,v) = - \a\ 2 - (ki + k )<0. 

JMouM! Jc* 

If K\ = cri(ui, v\) < then the Gauss curvature of Mq along C is negative. 
By Lemma 4.6 Mq is an unduloid or a catenoid. As 9{ < 7r/6 the vectors 1/1, 
which are tangent to the generating curve T± of Mq at their endpoints, are 
either horizontal or upper pointing. Therefore Mq contains a nodal region of 
the Gauss curvature in its interior, which implies that Mq is unstable [RR, 
Thm. 3], [PR, Prop. 4.1]. 

So for any component of R\ we have an admissible function u such that 
Q(u,u) < with support inside the boundary of the component. If we had 
three connected components in R± then we could get an admissible function 
satisfying the mean value zero property (3.11), which gives instability of the 
considered double bubble, a contradiction. □ 
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Remark 6.6. 

(a) An alternative, computational proof of Proposition 6.5, using Proposi- 
tion 6.1, is outlined by Heilmann et al. [HLRS, Prop. 4.6] (see [M2, 14.11, 
14.13]). 

(b) Reichardt et al. [RHLS] extended the arguments of Section 5 to prove 
the double bubble conjecture assuming only that one region is connected, 
thus providing an alternative to proving Proposition 6.5. 



7. Proof of the double bubble conjecture 

Theorem 7.1. The standard double bubble in M 3 is the unique area- 
minimizing double bubble for prescribed volumes. 

Proof. Let E be an area-minimizing double bubble. By Propositions 6.2 
and 6.5, either both regions are connected, or the region of larger volume 
and smaller pressure is connected and the one of smaller volume and larger 
pressure has two components. By Proposition 4.1, £ is either the standard 
double bubble or a bubble like the ones in Figures 8 or 9. As £ is stable, by 
Corollary 5.3 and Proposition 5.8 we conclude that it must be the standard 
double bubble. □ 

7.1. Immiscible fluid clusters. The methods of this paper extend to double 
clusters in which the three interfaces carry different costs, so-called immiscible 
fluid clusters (see M2, Chapt. 16]). We assume that the costs 001,002,012 
satisfy strict triangle inequalities, such as 

£02 = ooi + 0,12 ~ ^02 > 0. 

Theorem 7.2. For nearly unit costs, if the smaller region has at least 
37% of the volume, then the standard double cluster minimizes energy. 

Proof sketch. Proposition 6.1 has the following generalization to least 
energy: 

(7.1) 2E(vi,V2) > \- 1/3 e iA(v 2 ) + £02^1) + £12^(^1 + v 2 ). 

When the costs aij and hence the Sij are all 1, (7.1) reduces to Proposition 6.1. 
When they are near 1, the proof of Proposition 6.2 still shows that both regions 
are connected. Now the simple plane geometry of Corollary 5.3 shows that the 
cluster must be standard. □ 
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Remark 7.3. For general costs, even for equal volumes, it remains an 
open question whether the standard double immiscible fluid cluster minimizes 
energy. The above proof applies whenever we know that both regions are con- 
nected. (The more complicated plane geometry of Proposition 5.8 (or [RHLS]), 
for the case when the larger region is connected but the smaller region has two 
(or more) components, does not immediately generalize, because the generat- 
ing curves no longer meet at 120 degrees.) Unfortunately, for general costs, 
even for equal volumes, (7.1) does not imply both regions connected. 

We thank undergraduate Ken Dennison for raising this question. 
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